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ABSTRACT 



We present simultaneous XMM-Newton and INTEGRAL observations of the 
luminous black hole transient and relativistic jet source GX 339-4. GX 339- 
4 started an outburst on November of 2006 and our observations were under- 
taken from January to March of 2007. We triggered five INTEGRAL and three 
XMM-Newton target of Opportunity observations within this period. Our data 
cover different spectral states, namely Hard Intermediate, Soft Intermediate and 
High/Soft. We performed spectral analysis to the data with both phenomeno- 
logical and more physical models and find that a non-thermal component seems 
to be required by the data in all the observations. We find a hardening of the 
spectrum in the third observation coincident with appearance of a broad and 
skewed Fe line. In all spectral states joint XMM/EPIC-pn, JEM-X, ISGRI and 
SPI data were fit with the hybrid thermal/non-thermal Comptonization model 
(EQPAIR). While this model accounts very well for the high/energy emission ob- 
served, it has several drawbacks in the description of the lower energy channels. 
Our results imply evolution in the coronal properties, the most important one 
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being the transition from a compact corona in the first observation to the disap- 
pearance of coronal material in the second and re-appearance in the third. This 
fact, accompanied by the plasma ejection events detected in radio on February 
4 to 18, suggest that the ejected medium is the coronal material responsible for 
the hard X-ray emission. 

Subject headings: Black hole physics - stars: binaries (GX 339-4) - gamma rays: 
observations - accretion, accretion disks - radiation mechanisms: non-thermal - 
radiation mechanisms: thermal 



1. Introduction 



GX 339-4 (also called V821 Ara) is a recurrent black hole transient X-ray binary (i.e. 
this source und ergoes frequent out b ursts followed by ve ry faint states) with few episodes in 



quiescence (e.g. iKong et al.l (120001 )). iHynes et al.l (120031 ) determined, on the basis of optical 
spectroscopic observations, that this source has an orbital period of 1.7 days and a mass 
function of /(M) = S.SiO.SM© (since /(M)ocMsin^(i), note that this represents a lower 
limit to the b lack hole ma s s). T he distance and the spectral type of the secondary star 
are uncertain. iHynes et al.l (120041 ) established the distance to be 6-15 kpc depending on the 
spectral type of the secondary (a subgiant), being later than F spectral type. Therefore, 
this black hole transient is a low- mass X-ray binary (L MXB). Slightly later, on the basis of 
binary paramet ers determina t ion. IZdziarski et al.l (120041 ) suggested the distance to be >7kpc, 
consistent with iHynes et al.l (120041 ) . The inclination of this bi nary system is k nown to be 
low, from the small emission-line velocity amplitudes observed. IWu et al.l (120011) determine d 
a plausible value for the inclination of i = 15° from spectroscopic data. iMiller et al.l (l2004al ). 
Miller et al.l (l2004bl ) and iMiller I (l2006al ) found that fits of the relativistic Fe Ka line imply a 
very low value of the inclination (z = 20 — 30°) as well. The latter corresponds to the inner 
disk inclination, i.e. not the inclination of the binary plane. 

The timing and spectral properties of GX 339-4 required a re-interpretation of the stan- 
dard classification schem e of the black hole s tates. The behavio r of G X 339-4 supports the 
view that, as reviewed by lHoman et al.l (120011 ) and Ivan der Klis I (120061 ) , X-ray states are irn - 
portant, but are not a simple function of the luminosity. As claimed by lHoman et al.l (120051 ). 
another variable apart from the mass accretion rate, M, drives black hole state transitions 
and at least two factors m ust drive state trans itions; the second factor may be, e.g., the 
compa ctness of the corona (iHoman et al.l (l200ll )) . Black hole X-ray transients an d related 
states jTanaka et all Jl995h : lchen et al. I Jl997h : iMcChntock fc Remillardl J2006h ) provide 



an excellent opportunity to investigate the importance of irradiation and non-thermal pro- 
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cesses in the X-7-ray spectra of X-ray binaries. Black hole binary outbursts typically begin 
and end in the low/hard state, with a high-energy spectrum that can be roughly described 
by a powerlaw with a spectral index F varying in the range 1.4-2.1. At some point in the 
outburst (coinciding with the presence of radio outflows) the source transits from the hard 
intermediate (HIMS) to soft intermediate (SIMS) states, characterized by a quenching of the 
radio emission (i.e. the so-called jet line). Then, the thermal di sk component lu r ninosit y 
increases and the h ig h energy powerlaw is soft er (T = 1.5 — 2.5) ( Miyamoto et al.l (Il991al ). 



Belloni et al.l (119971 ) , iHoman fc Belloni I (120051 ) , see Ivan der Klis I (120061 ) for a review) . The 
transition between th ese states occurs very rapidly; they were found in EXOSAT observa- 
tions of GX 339-4 bv lMendez fc van der Klisl Jl997h. Th e Very/High (hereafter VH) state 
was first observed in GX 339-4 (IMiyamoto et al.lll991al ). when very peculiar X-ray vari- 
ability was dete c ted wit h Ginga/LAC (6 Hz QPO, 10-60s duration dips and "fiip-fiops"). 
Miyamoto et al.l (Il991bl ) even proposed a jet model in order to account for the observed 
time variability in this state, on the basis of the s ize of the Comptonizing coronae inferred 



from observations of GX 339-4 in the VH state (Mivamoto et al.l 1991a). Recent studies 



(iHoman et all (120011 ). iMendez fc van der Klis I (Il997r ). iKlein-Wolt fc van der Khs I (120071 )) 
point out that the VH state is really an intermediate state but at a different luminosity 
level. 



Since its diskovery more than 30 years ago by lMarkert et al.l (119731 ) . the source has been 
extensively studied with several optical, infrared. X-ray and 7-ray observatories. At higher 



energy its emission has 



Grebenev et al 



1993|) 



)een studied by all th e major observat ories, i.e., GRAN AT / SIGMA 



Bouchet et al 



(Il993h ). CG RO/OSSE iG^abelskY eTaDll995l), CGRO/B ATSE 



1994[) iRubin et al.N l998hl. ASCA Jwilms et al.Nl999h). (7moa JEbisawa et al. 



Harmon et al. , _ . .. . . ^ , 

Jl993h and lUeda et all Jl994h l RXT E (e.g. Izdziarski et al.l J2004h. 



Bel 



and INTEGRAL Jjoinet et all J2007h . iBelloni et all J2006h l iHoman et al 



oni et al. 



(119991)) 



(I2OO5I ) presented 



a multi-wavelength study of the 2002 outburst and suggested a non-thermal jet origin for 
the optical/near-infrared emission in the low/hard state, while the accretion disk dominates 
in the high/soft state. This source is very intriguing in the optical, since its photomet- 
ric optical behavior shows a very low amplitude modulation superposed by some flickering 
events. In addition to the QPOs detected in X-rays (usually appearing in the low/hard 
state and maintained in the intermediate states, while quenched in the high/s oft state; see 



(119821 ) 



Be 



loni et al.l (120051)), this s ource also shows optical QPOs (f ^0.064 Hz; iMotch et al. 



Steiman- Cameron et al.l (119971 ) ) . A very interesting issue is also the occurrence of 



optical states corresponding to the X-ray states (IMotch et al.l (Il983l ). iMotch et al.l (119851 )). 
This kind of variability is difficult to explain if optical radiation comes from reprocessing of 
X-rays and seems to demands optical activity originated in the base of a jet. 



GX 339-4 may harbor a black hole with a high spin parameter. iMiller et al.l (l2004al ) 
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and iReis et al.l (120081 ) found very small inn er di sk radius of 2 — 3 Rg (where Rg = GM / 
and a = cJ/GM'^; see iBardeen et al.l (Il972[) an d iThorne I (Il974j )). implying a dimensionless 
spin parameter of a > 0.8 — 0.9 jMiller" 2006a). Re lativistic radio jets with v/c > 0.9 have 



recently been observed, detected (iGallo et al. 



20041 ) during a transition from the low/hard 



to the high/soft state, therefore increasing the list of known microquasars showing similar 
radio and X- ray properties . A rad io jet is usually reported in the low/hard state of GX 339- 
4, for which iFender et al.l (119991 ) reported a quenching of the radio emission during the 
high/soft state (by a factor of >25 compared to the l ow/hard stat e ). Th e low /hard stat e 
would be associate d by a steady state of the outflow (ICorbel et al.l (120001 ). iFender I (120011 ) . 



Gallo et al. I (120031 )). It was also proposed that during state transitions the radio eniission 
results in one or more discrete ejection events, as already observed by iGallo et al.l (120041 ) 
during their detection of jets. GX 339-4 shares X-ray timing and spectral properties with 
the classical black hole Cyg X-1, although the former exhibits more f requent state changes 



and a l arger dynami c rang e of soft X-ray luminosity (ITanaka et al.l (119951 ) iZdziarski et al. 
Jl998h . iNowak et all Jl999h . 



While the physical nature of the soft component is commonly associated with an op- 
tically thick accretion disk, there is no consensus about the origin of the hard powerlaw 
component, where a powerlaw is probably a simplification of a much more complex reality. 
The energy spectra include also additional components, important fo r the physics of accre- 



tion, s uch as the Fe Ka fluorescence feature (see, e.g., iMiller I (120071 ) . iReynolds &: Nowak 



(120031 )) and Compton reflection bump, both being different aspects of the same physical ori- 
gin (i.e. fluorescence of the Fe ions and Compton back-scattering, respe ctively, both being 



the m ost obvious reactions of an irradiated disk by a high-energy source; iGeorge fc Fabian 



(Il99ll )). An extremely high ionisation may inhibit the formation of emission and absorption 
lines in the inner disk, which would explain the lack of Fe fluorescence emission line 
in some systems. One important parameter to measure is the presence/absence of a high 
energy cut-off in the spectrum, for which observations at energies > 20k e y are necessary. 
It is known since a long time (jSunyaev fc Titarchuk I (ll980l ) ]Grove et al.l (Il998l )) that the 
spectrum in the low/hard state show a cut-off around 100 keV. This cut-off was interpreted 
as being a consequence of the thermal (inverse) Comptonization of photons from the disk on 
a distribution of electrons thought to represent a corona (for which the geometry is a matter 
of debate). 

In this paper we report on spectral analysis of simultaneous INTEGRAL and XMM/EPIC- 
pn data, during five and three Target of Opportunity observations of GX 339-4, respectively, 
occurring between 2007 January 30 and March 31 (epochs 1-5 below; see Table [T] for more 
details). We fit the obtained spectra with phenomenological models in Section l3?n and we 
observed that the source spectral states evolved from one observation to the other, as re- 



- 5 - 



ported in Section [3^ Fits of the spectrum with both phenomenological and Comptonization 
models are presented in Section I3.H [5] and H] and we end with the discussion of the resuhs 
in Section O 



Observations 



In 2006 November, X-ray activity of GX 339-4 was detected with the Rossi X-ray 
Timing Explorer (RXTE; ISwank et al.l (l2006l )). The source had an almost constant flux 
until end of December 2006, when the hard (15-50 keV) X-ray flux increased by a large 
a mount. It reached its brightest level since 2004 November, as detected by SWIFT/BAT 



( iKrimm et al.ll2006l ). At the end of J anuary 2007, observations were triggered with the 
INTEGRAL satellite (IMiller et al.ll2007l ). During the following months, the source underwent 
an evolution from the low/hard to softer states, as shown in preliminary results reported in 
Caballero-Garcia et al. 2007a-e. Observations in radio with ATCA during February 4 to 



18 re vealed that the source was undergoing a series of plasma ejection events ( ICorbel et al. 



20071 ). Our INTEGRAL observations cover 5 Target of Opportunity observations of f5il30ks 
each, spread from 2007 January 30 to March 31 (see Tabled]). We also carried out a series 
of 3 ToO observations with the XMM-Newton satellite, simultaneous with the second, third 
and fifth INTEGRAL observations. Details about the dates and exposure times appear in 
Table [D 



2.1. INTEGRAL 



The data described were obtained with INTEGRAL, usin g the following instruments: 
the spectrometer on INTEG RAL fSPI: IVedrenne et all jSosI)), the INTEGRAL IBIS Soft 
Gamma- Ray Imager ( ISGR I; iLebrun et al.l (120031 )) and the Joint European X-ray Monitor 
(JEM-X; llund et al.l J2003h . SPI, ISGRI and JEM-X allow images to be obtained through 
the coded mask technique. ISGRI is optimized for 15 keV to 1 MeV imaging and SPI is 
optimized for high-resolution spectroscopy in the 18 keV to 8 MeV band. The former provides 
an angular resolution of 12' full-width half maximum (FWHM) and an energy resolution, 
E/AE of ^12 (FWHM) at 100 keV. SPI provides an angular resolution of 2.8° (FWHM) 
and an E/AE of 430 FWHM at 1.3 MeV. JEM-X has a fully coded Field of View (FOV) of 
4.8° diameter and an angular resolution of 3' FWHM, with a spectral resolution of 1.2 keV 
(FWHM) at 10 keV. JEM-X has a spectral resolution of E/AE = 0.2 at 6 keV, thus providing 
medium resolution spectral capabilities in the energy range of 3-35 keV. 
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The difference in the (effective) exposure times between the INTEGRAL instruments 
given in Table [1] (JEM-X, ISGRI and SPI) are due to the difference in the dead times and 
variation in the efficiency along the fields of view. The dithering pattern used during the 
observations was 5x5 (square of 25 pointings separated by 2.17 degrees centered on the 
main target of the observation); this is the best pattern in order to minimize background 
effects for the SPI and ISGRI instruments in crowded fields. Data reduction was performed 
using the standard Off-line Science Analysis (OSA) 7.0 software package available from the 
INTEGRAL Science Data Centre (ISDC). 

Because of the steep fall in response off-axis in the case of JEM-X, and because of its 
reduced FOV (~5° of diameter), we limited the pointing radius with respect to the GX 339- 
4 position to be within 4°. In the case of SPI and ISGRI, with large Fully Coded Fields 
of Views (FCFOV) (16°xl6° for SPI and 9°x9° for ISGRI), pointing selections were not 
necessary. In total, 149, 254 and 107 pointings (or Science Window - each having exposure 
times lasting from 1800 to 3600 s) data were used for SPI, ISGRI and JEM-X, respectively. 
The difference in the number of pointings between SPI and ISGRI is due to the fact that 
in epoch 4 and 5 the source was too weak in hard X-rays {E > 20keV) to be detected 
with SPI. Due to the fact that SPI is a high resolution spectrograph in its energy range 
(20-8 000 keV), it is not optimized for the detection of sources without taking into account 
previous information about the spatial distribution of them. Thus, we used the list of found 
sources obtained with ISGRI in the (20-40) keV mosaic images (see Figure [1]) as our input 
catalog. This resulted in GX 339-4 as the only source detected in the 200-300 keV energy 
range (Figure [2]) with a detection significance of 4. 80". 

We combined single revolution SPI spectra from epochs 1-3, considering that there was 
no significant evolution within each revolution (as shown by the relative constancy of the 
INTEGRAL and XMM-Newton light curves, except for epoch 3 -see below-; Figures [3] and 
H]), in order to increase the signal to noise ratio. We obtained 3 SPI spectra in total. The same 
was done for the low energy instruments (JEM-X and ISGRI), thus obtaining five spectra, 
namely one for each epoch (see Tabled]). We applied 2% systematic errors to the JEM-X, 
ISGRI and SPI spectra. We restricted our analysis in the energy ranges of 5 — 20, 18 — 200 
and 23 — 300 keV for the JEM-X, ISGRI and SPI spectral analysis, as recommendecQ. The 
SPI and ISGRI spectra were rebinned at high energies (>200keV) with the FTOOL grppha 
procedure to reach the detection level of 2a per spectral bin. 



^For OSA software download, cookbook documentation and cross-calibration issues, we refer to the web 
page: |http://isdc.unige.ch/?Soft+download 
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2.2. XMM/EPIC-pn 



The XMM-Newton Observatory (jjansen et al.l 1200 ll ) includes three 1 500 cm^ X-ray 
telescopes each with an European Photon Imaging Camera (EPIC; 0.1-15 keV ) at the focus. 
Two of the EPIC imagiiig spe ctrometers use MOS CCDs (ITurner et al.ll200ll ) and one uses 
pn CCDs (IStriider et al.ll200ll ). The Reflection Grating Spectrometer, (RGS; 0.35-2.5 keV, 
Den Herder et al.ll200ll ) are located behind two of the telesco pes. In addition, th ere is a co- 
aligned 30 cm diameter Optical/UV Monitor telescope (OM. iMason et al.ll200ll ). providing 
simultaneous coverage with the X-ray instruments. 

GX 339-4 was observed by XMM-Newton for 16.8 ks on 2007 February 19 (hereafter 
called Obs 1 and coinciding with INTEGRAL Epoch 2), for 17.8 ks on 2007 March 5 (here- 
after called Obs 2 and coinciding with INTEGRAL Epoch 3) and for 20 ks on 2007 March 
30 (hereafter called Obs 3 and coinciding with INTEGRAL Epoch 5). We refer to Tabled] 
for details. 

The thin optical blocking filter was used with the EPIC camera. The EPIC pn camera 
was operated in burst mode due to the high count rate of the source (from >6700 s^^ to 
>3000 s~^ for Obs 1 to 3, respectively). For the same reason, the EPIC MOS cameras were 
not operated. The individual RGSl and RGS2 CCD chips were read out sequentially. This 
reduces the frame time from 4.8 to 0.6 seconds (4.8/8) and consequently the pile-up by almost 
an order of magnitude. All the X-ray data products were obtained from the XMM-Newton 
public archive and reduced using the Science Analysis Software (SAS) version 7.1.0. In pn 
burst mode, only one CCD chip (corresponding to a field of view of 13^6x4^4) is used and 
the data from that chip are collapsed into a one-dimensional row (4^4) to be read out at high 
speed, the second dimension being replaced by timing information. The duty cycle is only 
3%. This allows a time resolution of 7 /is, and photon pile-up occurs only for count rates 
above 60 000 s~^. Only single and double events (patterns to 4) were selected. Ancillary 
response files were generated using the SAS task arf gen. Response matrices were generated 
using the SAS task rmf gen. In this paper we do not use RGS spectra because our main aim 
is to make a study of the continuum properties of the spectra. 

PN Charge Transfer Inneficiency (CTI) rate dependency has been observed for count 
rates > 200 — 300cts/s, i.e. for very bright sources CTI correction overpredicts the CTI 
losses. This can result in an up to 2% apparent gain shift most visible for spectral features 
like edges or lines §. Hereafter, we applied 2% systematic errors to the EPIC-pn spectra 



^For details see: http://xmin2.esac.esa.int/docs/documents/CAL-TN-0018.pdf (EPIC Calibration Status 
Document), page 4. 
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reported below. We used EPIC pn data in the 0.7-10 keV energy range. 



2.3. Light curves 

In Figure El we show the SWIFT/BAT and the RXTE/ASM light curve in the 15- 
50keV and 2-12 keV energy ranges during the entire outburst. The horizontal lines shown 
in the Figure indicate the time intervals (150 ks each, in five observations) over which our 
average spectra were obtained. Not e that black ho l e transients u s ually b egin and end their 



outbu rsts in the low/hard state (see iNowak I (119951 ): iFender et al.l (120041 ): iHoman fc Belloni 



(120051 )). and the 2007 outburst of GX 339-4 is no exception. As can be seen in Figure El the 
beginning of the outburst started earlier at high energies (15-50 keV) rather than in softer 
X-rays (2-12 keV), when observations of epoch 1 were made. The other four epochs were 
made during a period when the source spectrum was much softer. It is interesting to notice 
the sudden increase of 15% fiux (with respect to the maximum value) of the 15-50 keV light 
curve, during the transition from second to third epoch observations, possibly associated 
with changes in the high-energy emission source, as we will discuss in Section |6l 

In Figure [3] we show the XMM/EPIC-pn light curves in the 0.6-10 keV energy range 
and calculated hardness ratio (counts in 2-10 keV band divided by counts in 0.6-2 keV) with 
a time binning of 120 s. We do not see spectral evolution during the time of each observation 
due to the constancy of both the light curves and hardness ratio (except for a small increase 
in the fiux during a short time period in epoch 3, but of the order of <6%). 



In Figure H] we show the ISGRI (20-200 keV) light curves during all the observations. 
The source was significantly detected during all the epochs with decreasing fiux (except for 
epoch 3, with higher fiux than for epoch 2). Within every epoch the fiux was relatively 
constant, except for epoch 3, in which the fiux was decreasing with time, with a total fiux 
variation of 70% around the mean fiux value. In spite of this, we combined all the revolution 
data to get the spectrum in this epoch as well (as did for the remainder epochs) in order to 
get higher signal to noise ratio. 



3. Analysis of spectra 



We performed fits to the combined JEM-X, ISGRI and, when avai lable, SPI and X MM/EPIC- 
pn spectra, for each of the five epochs (see Tabled]) using XSPEC ( lArnaud Ill996l ) v. 11. 3. 
All errors quoted in this work are 90% confidence errors, obtained by allowing all vari- 
able parameters to fioat during the error scan. The hydrogen column density was set to 
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be free when XMM/EPIC -pn data were available and fixed at Nh = 0.5x10^^ atom cm ^ 
( iDickey fc Lockman 1119901 ) when only INTEGRAL data were available. 

Our main aim in these fits is to characterize the broad continuum thanks to the broad 
XMM-Newton and INTEGRAL spectral coverage. To account for uncertainties in relative 
instrument calibrations, we fixed JEM-X multiplicative calibration constant to 1 and let that 
of ISGRI and SPI free to vary in the fits for the different data sets as shown in Table O 
In the case of using both INTEGRAL and XMM/EPIC-pn spectra, we fixed the EPIC-pn 
multiplicative constant to 1 and let that of JEM-X, ISGRI and SPI free to vary, in order 
to account for cross-calibration uncertainties as well|f|. The presence of an XMM/EPIC-pn 
instrumental line was clear in the spectra of epochs 2 and 3 and we added a gaussian emis- 
sion line centered at 2.28 keV and 2.25±0.03keV, respectively. For epoch 5 two absorption 
instrumental lines appeared, centered at 1.82 keV and 2.14 keV, respectiv ely. As mentioned 
in the XMM-Newton analysis of GRO J1655-40 JPiaz Trigo et al.lboOTl ): where these fea- 
tures were present as well, structured residuals near 2.3 keV are probably due to an incorrect 
instrumental modeling of the Au mirror edges, whilst those near 1.8 keV are probably due 
to an incorrect modeling of the Si absorption in the CCD detectors. 



3.1. Fits of joint XMM/EPIC-pn and INTEGRAL spectra with 

phenomenological models 

We performed preliminary fits to the spectra of every single epoch with t he phenomeno 



l ogica l model consisting of a multicolor disk black body (diskbb in XSPEC) of iMitsuda et al. 



( 1l984j ) plus a powerlaw convolved by absorption at low energy in order to describe the high- 
energy emission. In none of the observations this model provided a good fit to the data 
(xl = 46,3.05,4.2,1.8,3.3 for u = 96,150,154,46,108 for epochs 1, 2, 3, 4 and 5, respec- 
tively). 

In epoch 1 the situation improved dramatically {xl = 2.76, u = 94) when a high-energy 
cut-off was included in the model (i.e. constant*phabs (diskbb + powerlaw) highecut; 
Ec~70keV). Nevertheless, some positive and negative residuals still remained at 6-10 keV 
and 10-20 keV, in the form of a steepening of the spectrum and of an edge, respectively, 
indicating the likely presence of a broad FcKq, line. We included a broad Fe gaussian line 
centered at Epo = 6.4 — 6.97keV, getting the final statistics of xl = 1-4 {i^ = 91). With the 



■^We restricted the values of the multiphcative constants to be 1. Ril±20%, 1±30%, 1±60% for 
XMM/EPIC-pn,JEM-X, ISGRI and SPI relative to XMM/EPIC-pn in order to obtain rehable fits. When 
XMM/EPIC-pn data was not available, the same ratio between the constants was preserved. 
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model constant*phabs(diskbb+gaussian+powerlaw)highecut, the final values obtained 
for the photon index and the folding energy are of F = 1.46j| n 'ni^ an d E = 66.4:~^^'^keV, 



respectively. These values resemble those found by iBelloni et al.l (120061 ) during the HIMS of 
the 2004 outburst. 

In epoch 2 low-energy positive residuals in the form of an excess near 1 keV were detected 
in the XMM/EPIC-pn spectrum and the fit greatly improved when considered {x1 = 1.9, 
u = 144). The obtained equivalent width is Wikev = 11.6±2.1eV. This feature is detected 
in a number of X-ray binaries and has been pr eviously modele d eithe r as an emissi o n line , 



or as an edge, and its nature is uncl e ar (e .g. iKuulkers et al.l (119971 ): ISidoli et al.l (120011 ): 



Boirin et al.l (120031 ) . iDfaz Trigo et al.l (120071 )). However, since this has not been detected 
before in GX 339-4 we think that this line probably is an instrumental effect due to incorrect 
modeling of Ne IX and Ne X and Fe L in the CCD detector. 

In epoch 3 positive and negative residuals appeared, centered at ~6 — 10 keV and ~10 — 
20keV, respectively (see Figure [6]). These residuals are compatible with a broad and skewed 
FeKc line and the respective Fe edge. As in epoch 1, we added a gaussian line centered 
at Epe = 6.4 — 6.97 keV and the fit improved {xl = 3.3, u = 146), but some residuals in 
the red wing of the Fe line (3-5 keV) still remained. These residuals point out to the Fe Kq, 
line having a relativistic origin (i.e. originate d very close to the black hole). Thus, we 
substituted the gaussian component by a laor (iLaor Ill99ll ) model, which includes general 
relativity effects in the modeling of the X-ray emission in the inner boundaries of a rotating 
black hole. The parameters of this model are the line energy E, the disk emissivity index 
q (where the emissivity is assumed to have a power-law form of J{r)(xr~'^; where g = 3 is 
expected for a standard disk), the inner radius of the line emission region in units of GM/c^ 
(with a lower limit of 1.235 Rg for a maximally rotating black hole), the outer line emission 
radius in units of GM/c^ (fixed to the fiducial value of 100 i?^), the disk inclination, and 
the line normalization. With this model, we obtained significant better statistics {xf, = 2.7, 
u = 143) and flat residuals in the 3-5 keV energy range, where the relativistic effects of the 
Fe complex are expected to appear 0. The inner radius of the disk, the emissivity index, 
the disk inclination and the normalization were free to vary in the fit. We allowed the value 
of the inclination to vary in the range 10°-60°. The value obtained for the line energy is 
very high (6.83±0.14keV), indicating a very ionized medium. The obtained value for the 
emissivity index is rather steep (g = 5.2l5;o)) being close to the value found in the VH 



"'Whilst in epoch 1 there are residuals compatible with the broad Fe line, we can not test if its origin 
is relativistic, due to the JEM-X lower energy limit of 5 keV, where the relativistic red wing is expected to 
appear. Thus, complex models dealing with the sophistication of the red wing line (as laor), can not be 
applied to the data of epoch 1. 
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state (g = 4.8 — 6 by iMiller et al.l (l2004a( )). The obtained value for the disk inchnation is 
i = 35'tj°- High values for the emissivity index indicate that the primary source is very 
close to the black hol e. Steep emis s ivity in dices (~7) have bee n observed from a number of 
Seyfert galaxies (e.g. 
was rather low (r^ = 



Fabian et al.l (120041 ): iFabian et al.l (120051 )). The inner radius obtained 



2.11±0.11 i?g), thus being clearly different from 6Rg ('exp ected value 



for a non-rotating Schwarzschild black hole) and close to the value found by iReis et al. 
( I2OO8I ) (Tin = 2 — 2.1 Rg). The presence of negative residuals at lOkeV led us to consider the 



effects of a smoothed absorption of the Fe ions in the inner disk. We included a smoothed 
absorption edge (smedge in XSPEC) and the fit improved, getting the final statistics of 
xl = 1-3 {u = 141). Relativistic smearing is expected to appear in the inner parts of the 
accreting disk and, with Fe line, accounts phenomenologically for reflection. This resulted 
in an edge energy of -Egmedge = 8.53±0.20 keV, compatible with He-like configuration for Fe 
and r = l.OQj^n 1^1 w i th a s meared edge width frozen to W^smedge = 5keV (close to the value 
found bv lReis et aP J2008h in the VH state; i.e. W^medge = 4.5^2:0 keV). 



In epoch 4 some residuals remained in the form of a smeared edge in the 10-20 keV en- 
ergy range. When considered (i.e. model constant*smedge*phabs(diskbb + powerlaw)), 
the fit improved {xl = 1-2, u = 44), meaning an F-test probability of 1.337E-04 that the 
improvement was made by chance. The value obtained for the threshold energy and the 
maximum absorption factor at threshold are i?smedge = 9.53lo;gkeV and r = 1.6±0.5, re- 
spectively, compatible with H-like configuration for Fe. 



Finally, in epoch 5, residuals in the form of an excess at >50keV and a smeared edge 
at 7-10 keV {xl = 3.3, u = 108) were present. In the residuals there was an emission 
line like in epoch 2 and we included an emission gaussian line centered at ~1 keV {xl = 2.4, 
z/ = 105). The equivalent width obtained was Wikev = 29±5 eV, thus higher than in epoch 2. 
When adding the smoothed edge (i.e. model constant*smedge*phabs (diskbb + powerlaw 
+ gaussian)) the fit greatly improved {xl = 1-3, u = 104). The obtained values for both the 
threshold energy and maximum absorption factor at threshold are -Egmedge = 9.22±0.25keV 
and r = 3.0±0.5, respectively, the former compatible with H-like configuration for Fe i The 
value for the absorption factor at threshold is very high and would indicate that it is not 
due to the presence of an absorber, rather being due to reflection effects. The excess above 
50 keV would indicate an important source of Comptonization by non-thermal particles, as 
occurring in epoch 1, as shown below. 

The values obtained for the best fit parameters with the best phenomenological fit are 
shown in Table [2] and the spectra with the best fit results are shown in Figures [7] and [HI 



'"^The edge for neutral iron is at 7.1 keV, for He-like iron at 8.8 keV, and for H-like iron at 9.3 keV. 
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In Table [3] we summarize the values for the photon indices of the powerlaw (analyzed in 
this Section) and un-absorbed luminosities of both the powerlaw and the disk (assuming a 
distance to the source of 8.5 kpc). 



3.2. State classification and related behavior 



The main features of the spectrum of epoch 1 are a rather flat photon index (F = 
1.461^002) and the clear presence of a high-energy cut-off at 66. 4j]?'? keV. The value found for 



the cut-off is very close to the value reported by iBelloni et al.l (120061 ) {Ec = 70lil^ keV) dur- 



ing HIMS, although the value found for the photon index in their case was clearly steeper 
(r = 1.92±0. 5). We conclude th a t the s pectrum of epoch 1 corre sponds in luminosity 
teelloni et aD Jl997h . IPender et aP kooi \. iHoman fc Belloni I koO^ ) and spectral char- 
acteristics to the HIMS. During epochs 2 and 3 the spectra are characterized by a steep 
high-energy emission ( rR::!2.33±0.02, 2.659 j ;n m'^ for epochs 2 and 3, respectively). In the 



classification scheme of iHoman fc Belloni I (120051 ). both correspond to the SIMS. The value 
of the photon index fr om epoch 2 re s emble s those reported previously by lNespoli et al.l (120031 ) 
(r = 2.44±0.10) and lBelloni et all J2006h (F = 2.3) during SIMS of GX 339-4 as well. The 



inner disk tempe ratures found are 0.8 91±0. 003 keV and 0.852^ynn« k eV, thus resembling the 



values found bv iMalzac et al.l (l2006f ) and ICadoUe Bel et al.l (120051 ) (i.e. 1.2±0.2keV) for 
Cyg X- 1 in Intermediate st ates as well. Observations of GX 339-4 during the 2004 out- 
burst by lMiller et al.l (l2004al ) during the VH state show the same photon index value for the 
high-energy emission (i.e. T = 2.611q;o^) as us for epoch 3, accompanied by a detection of a 



relati vistic Fe line, as well. Indeed, in the most recent state prescription (IHoman fc Belloni 
20051 ). the VH state is the same as the SIMS. 



The spectrum during epochs 4 and 5 is characterized by a steep high-energy emission 
(r = 2.39±0.06, 2.65±0.02). As can be seen in Figure [3] the hardness ratio (defined as 
HR=(2-10) keV)/(0.2-2) keV) of the light curves are very soft in epochs 4 and 5. The inner 
disk temperatures found for these epochs are 0.75lQ Qg keV and 0.698±0.003 k eV, respectively. 



thus r e sembling the va l ues fo und in the HS state for Cyg X-1 (i.e. O.GkeV iGierlihski et al. 



(1l999l ). iFrontera et al.l (120011 )). The softness, the inner disk temperatures, the values for 
the photon index and the high value for the disk normalization component suggest that 
the source in those epoch s is in the high/soft state, considering the classification scheme of 
Homan fc Belloni I (120051 ). The spectral evolution is presented in Figure [HI 



As explained above, our observations of GX 339-4 cover the spectral evolution from 
the HIMS to high/soft state, passing through the SIMS. As shown in Section [3.11 and in 
Table [31 this spectral evolution involves very significant spectral changes which are strongly 
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correlated. The changes in the spectral shape are strongly correlated with the ratio of the 
powerlaw to the disk luminosity. In the first epoch there was almost no disk contribution 
in the spectrum. It showed a cut-off at Ec = 66.4t??keV in the powerlaw. The cut-off 
disappeared in the second epoch and the powerlaw emission decreased notably, being just 
1/7 of the disk luminosity. In the third epoch the spectrum hardened again considerably (the 
powerlaw luminosity being 1.7 times higher than the disk luminosity, the latter decreased 
by 50% with respect to the second epoch). The total luminosity increased by 50% in the 
transition from epoch 1 to 2 and remained more or less unchanged during epochs 2 and 3, 
while decreasing in epoch 4 and 5 (again decreased by 50%). The luminosities obtained for 
all the epochs were within the range of 20 — 50% the Eddington luminosity (assuming the 
compact object being a black hole of 6M0). 



4. The high-energy excess 



During the LH state of the GX 339-4 outburst occur r ing in 



(above 200 keV) was observed with OSSE (jjohnson et al. I (119931) 



991, a high-energy exce ss 



Wardzihski et al.l (120021 )1 



In 2004, during the increasing phase of the 2004 outburst, iJoinet et al.l (120071 ) observed a 
similar feature with SPI when GX 339-4 was in its canonical LH state. Similar features 
have been recently reported for other b lack hole transi e nts in LH states observed d uring 
the decreasing phase of their outbursts (iKalemci et al. 1 (|2006[),l Kalemci et al. ] jSoi)). It 
is also more commonly assoc iated with state transitions ( Malzac et al. 20061 ) and with soft 
states (iGierlihski et al.lll999l ). 



Such a hard tail is ussually attributed to the presence of a small fraction of non-thermal 
electrons in the hot Comptonizing plasma. In the HIMS of our observations (epoch 1), GX 
339-4 shows a highly significant cut-off at 66.4l5'i keV. This indicates that soft photons are 
Comptonized by a thermal electron population with a temperature kTe that is related to 
the cut-off energy. Both for the HIMS (epoch 1) and SIMS (epochs 2 and 3) non-negligible 
values of inth/^h were found (as shown below in Section E]). 

In order to investigate further the possible presence of a non-thermal component dur- 
ing the HIMS and SIMS, IBI S and SPI spectra have been fit with the COMPPS model of 
Poutanen fc Svensson I (119961 ) (see Table Hj). In this thermal Compton model, the spectrum 
is computed for a homogeneus spherical hot corona, inside which blackbody photons are in- 
jected and then upscattered by hot Maxwellian electrons. The reflected spectrum is smeared 
out by rotation of the disk due t o special and gen eral relativistic effects using diskline-type 
kernel. As in the comptt model (ITitarchuk Ill994j ). the corona is parametrized by its electron 
temperature, Te, and the radial Thompson optical depth, r = ngaxH, where rif. is the elec- 
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tron density, ar is the Thompson cross section and H the radius of the sphere. Using this 
thermal Comptonization model, some residuals are present at high-energy, especially in the 
spectrum of epoch 1. In order to mimic the presence of a non/thermal component, we added 
a powerlaw to the pure Comptonization model (see Figure [TU]) . This led to an improvement 
of the that is significant for epoch 1 (see Table Hj). The F-test probability that this im- 
provement was made by chance is p = 6A61E — 02 for epoch 1. On the contrary, because 
of the low statistics in epochs 2 and 3, adding a powerlaw to the thermal Comptonization 
model is not required. However, we notice that the electron temperatures obtained with the 
thermal Compton plus powerlaw model are rather close to the values obtained with EQPAIR 
fits for all the epochs. 



5. Fits of joint XMM/EPIC-pn and INTEGRAL spectra with the hybrid 
thermal/non-thermal Comptonization EQPAIR model 

5.1. The model 

We saw in the previous Section that a high-energy excess is visible in the spectra. 
Since this high-energy excess can not be well described by thermal Comptonization, we 
go into more sophisticated models containing a hybrid distribution of thermal/non-thermal 
particles. From the variety of Comptonization models available in the l iterature, we have 



chosen the hybrid thermal/non-thermal Comptonization EQPAIR model (ICoppi Il2000l ). and 
relativistic line emission and disk reflection. The best-fit values for the parameters for all 
the epochs are shown in Table O Because the presence of an instrumental line at ~1 keV in 
the XMM/EPIC-pn spectra of epochs 2, 3 and 5 was analyzed in Section 13.11 this will not 
be discussed further. 

In the EQPAIR model, emission of the disk/corona system is modeled by a spherical 
hot plasma cloud with continuous acceleration of electrons illuminated by soft photons from 
the accretion disk. At high-energies the distribution of electrons is non-thermal, but at low 
energies a Maxwellian distribution with temperature kT^ is established. 

The disk spectrum incident in the plasma is modeled as coming from a pseudo-Newtonian 
accretion disk (diskpn in XSPEC) extending from Rout = 100 -Rg down to the minimum 
stable orbit Rout = QRg. Its spectral shape is then characterized by the maximum color 
temperature of the disk, kTpn. Previous observations of Cyg X-1 with X-ray telescopes in- 
dicate temperatures ra n ging from 0.1 keV in the LH state up to 0.6 keV in the HS state 



(IGierhhski et al.l (119991 ) . iFrontera et all fl200l[ )l 



The model (EQPAIR) takes into account angle dependence, Compton scattering (up 
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to multiple orders), photon pair production, pair annihilation, bremsstrahlung, as well as 
reflection from a cold disk. This model in c ludes ionized Compton reflection (as the pexriv 
model in XPSEC; iMagdziarz &: Zdziarski I (jl995l )). The reflection component is smeared out 
by rotation of the disk due to special and general relativistic effects using d iskline-type 



kernel to relativistic motion in the d isk (as in th e diskline model in XSPEC; iFabian et al. 



( 119891 ) )■ We used the LAOR model (ILaor Ill99ll ) to model the relativistic iron line emission, 
with the emissivity index (/?) free and tied to the opposite value of that of the EQPAIR. 
The disk reflection component is calculated for neutral material with standard abundances, 
unless otherwise noted. We asumed free inclination and inner disk radius in the range of 
10° — 60° and 2 — QRg. The outer disk radius was fixed to 100 -R^. 

The EQPAIR model provides the injection of a non-thermal electron distribution with 
Lorentz factors between and T^ax and a powerlaw spectral index Tinj- The system 
is characterized by the power (i.e. luminosity) Li supplied by its different components. 
We express each of them dimensionlessly as a compactness parameter, ii = LiCt / (RTrieC^) , 
where R is the characteristic dimension and ctt the Thompson cross-section of the plasma. 
Thus, is, ith, i-nth and ih = C-th + inth correspond to the power in a soft disk entering the 
plasma, thermal electron heating, electron acceleration and the total power supplied to the 
plasma. The total number of electrons (not including e+e~ pairs) is determined by tt, 
the corresponding Thompson optical depth, measured from the center to the surface of the 
scattering region. If we consider injection from pairs e"*"e~, then the total optical depth of 
the thermalized scattering electrons/pairs is expected to be tt^>t^. 

In general, the spectral shape is insensitive to the exact value of the compactness, but it 
depens strongly on the compactness ratios {dh/C-s and inthjih)- Thus, we froze to a fiducial 
v alue {If, = 10), as co mmonly reported for other sources with similar characteristics (e.g. 



Ibragimov et al.l (120051 )). In the fits reported below, we fit the data with a powerlaw dis- 
tributed injection of electrons and compare with the results obtained with a mono-energetic 
distributed injection of them as well. The former distribution is expected in the case of shock 
acceleration of particles, while the second could b e achieved in reconnection events that are 
expected to power the corona (iGaleev et al.lll979l ). 



5.2. Fits to the data 



For the first epoch of our observations, due to the low disk contribution, we froze the 
disk temperature to a fi ducial value of kTpn = 300 e V, as an intermediate value between LH 



and HS values found by iGierlihski et al.l (119991 ) and iFrontera et al.l (120011 ) for Cyg X-1. The 
relativistic Fe disk line emission was significant (with a chance probability of 4.823E-05). 
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The most relevant factor is the inchnation angle of the reflecting material, clearly favouring 
high values. We obtained worse residuals and statistics {xl = 1-8, u = 82) for the case 
of i = 30° with respect to the case of i = 50° {xl = 1-0, = 82). Since this value is 
in conflict with the value obtained in Section 13.11 and the values obtained from fitting the 
X-ray data in the literature, we will adopt i = 30° in the fits of the remainder epochs B 
High values for the ionization are preferred by the data, with an improvement of Axl = 0.5 
(for 82 d.o.f) with respect to the neutral case. The spectral cut-off already noticed in 
Section 13.11 can be modeled here by C omptonization of the rmal distributed particles and 



reflection from a cold disk. As shown in lPetrucci et al.l (1200 ll ). with an electron temperature 

-0.15 
-0.18' 



of /cTg = 27.5±1.2keV and an optical depth of r = 2.39tQ;}g, it is possible to reproduce 



the spectral cut-off at 66.415'^ keV. The spatial distribution of the Comptonizing particles 
is compact {(-h/^-s = 3.9;!^q;3) and mildly non-thermal {inth/^h = 0.40!!^q;o3). The high-energy 
distribution is well fit with a powerlaw injection of non-thermal electrons with relatively low 
photon index {Tinj = 1.80^g;o5) with Lorentz factors within the range 1.3-100 (as considered 
in the remainder epochs as well). This issue means the clear presence of non-thermal particles 
in the plasma. When fitting with a mono-energetic distribution of electrons with a Lorentz 
factor of r = 5, the fit slightly worsened {xl = 1.1, = 83 and worse description of the 
highest energy emission) with respect to the previous fit (i.e. xl = 1-0; = 82). Albeit not 
significant, the spectrum seems to be better described by a powerlaw injection of non-thermal 
electrons. Clearly, higher energy data are needed in order to disentangle these models. 

In the second epoch, XMM/EPIC-pn data allowed us to constrain the disk temperature 
{kTpn = S97~^l'leV). The relativistic Fe disk line emission was found to be zero (no improve- 
ment of the statistics at all, when considered). The data preferred a neutral reflecting disk 
(with Ax^ = 0.15 of improvement for 142 dof) and the final statistics is xt — ^-"^ — 142). 
The statistics is poor and the disk temperature very high, which would indicate that EQPAIR 
model is not adequate enough to describe the softest part of the spectrum of this epoch 0. 
Nevertheless, the high energy spectrum is well described by a purely non-thermal dominated 
distribution of particles {inth/^h = 0.9±0.1). However, the photon index of the powerlaw 
injection distribution of the Comptonizing electrons is very steep {Tinj = 2.87^Q;og), which 
indicates the presence of mildly relativistic electrons. Both the coronal compactness and 



^We repeat the same analysis by considering i 
change significantly with respect to the case i = 30° 



50° for epochs 2, 3, 4 and 5 and the results do not 



^As reported in iMalzac et al. ( 200G ) and 



Cadolle Bel et al.l (j2005f ). such high inner temperatures are 



not realistic for Cyg X-1. iMalzac et al.l (j2006[ ) included an additional warm Comptonization model to the 
EQPAIR model to describe the soft part of the spectrum during the Intermediate state finding very good 
agreement with the data. However, this kind of work is out of the scope of this paper. 
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opacity were found to be very low (4/4 = O.OS^o:"^ and Tp < 0.02) flNowak et al.ll2002h . 
Albeit the spectrum does not show significant relativistic smearing, the reflection covering 
factor was found to be significant {Q/2n = 0.401o;q4). We obtained the same quality fit 
when considering a mono-energetic distribution of electrons with a Lorentz factor of F = 5 
{xl = l.7,u = U3). 

In the third epoch the disk temperature was f ound to be fcTn„ = 704±2eV, thus re 



sembling the values found in soft states for Cyg X-1 (iGierlihski et al.l (Il999l ). iFrontera et al. 



( l200ll )). Albeit in the phenomenological fits reported in Section [XT] the presence of a rela- 
tivistic iron line was clear in the residuals, we notice the lack of its significance in this model 
(just a chance probability of 1.151E-04 compared with 1.040E-09 from the phenomenological 
fits). This issue can be understood by the fact that the pexriv model (which is incorporated 
in EQPAIR to account for reflection) is not a proper description of the reflection occurring 
in the inner parts of an accretion disk, since it does not include the Compton scattering 
(and consequent blurring in energy) that must occur. The statistics obtained is xl = 1-5 
(z/ = 138) and the high-energy spectrum is well described taking into account thermal and 
non-thermal distributed particles in the plasma, by the same amount, as occurring in epoch 
1. The photon index of the powerlaw distribution of non-thermal injected particles is very 
steep {Tinj = 3.42^q;q4), meaning a mildly relativistic population of non-thermal electrons. 
Both the coronal compac t ness a nd the opacity {^h/^s = 0.28^q;o'^ and Tp = 1A3~^q^I) are as 



found in iDel Santo et al.l (120081 ). The relativistic smearing in this spectrum is important, 
meaning that reflection probably occurs in the inner regions of the accretion disk. High 
values for the reflection covering factor were favored and we froze this value (f2/27r = 1). 
We considered a mono-energetic injection of electrons case as well (F = 5) and we obtained 
clearly worse agreement with the data {xl = 1.9, u = 139) and with a large discrepancy of 
the highest energy bins. In this epoch a powerlaw injection of electrons is clearly favored 
against the mono-energetic injection scenario. 

The results obtained for the fits of epoch 4 resemble those obtained for epoch 3. The 
maximum disk temperature was frozen to a fiducial value of 500 eV, for the same reason 
as well. When considering a powerlaw injection of non-thermal electrons, the fit statistics 
obtained is xl = 1-5 {i^ = 44). The corona is relatively compact, mildly non-thermal and 
relativistic smearing is important, thus giving a very high value of the reflection covering 
factor {Q/2'K = 1), as in epoch 3. Again a mono-energetic injection of electrons is a very 
poor description of the data {xl = 4.5; u = 45) with a large discrepancy of the highest 
energy bins. It is worthy to notice as well that in this case the iron line emission (modeled 
here as a LAOR component) was not significant in the data and even worsened the fit when 
considered (by Axl = 0.3, 0.2 for cases of i = 30° and i = 50°, respectively). 
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For epoch 5 we could determine the disk temperature [kTpn = 550±5eV). The relativis- 
tic hne was significant (with a chance probabihty of 1.029E-06) and the refiection covering 
factor very high {Q/2tt = 0.72j;Q;];Q). The fit statistics was = 1-6 = 99). As in epochs 
1, 3 and 4, mono-energetic injection of non-thermal electrons is not a good description of 
the data, due to the worse statistics {xl = 1-7, = 100) and discrepancy of the high-energy 
channels. The observed temperature of the disk is cooler (HS state) than in epoch 3 (SIMS 
state). This fact is accompanied by the increase of the compacticity for the latter and would 
indicate that the corona becomes more compact as the disk temperature increases. 

Spectra and models of the five different periods have been plotted in Figure [91 



6. Discussion and conclusions 



During our observations, we detect spectral evolution from the HIMS to the HS state, 
passing through SIMS. We performed phenomenological and more physical motivated fits to 
the data and found that during the HIMS the spectral cut-off at 66A~^i'^keV is due to the 
thermal kinetic distribution of electrons in the plasma. For the SIMS and HS states we do 
not find a significant cut-off in the data but applying analysis with Comptonization models 
(EQPAIR and COMPPS) we detect a high-energy tail in all the epochs that can be under- 
stood as evidence of a non-thermal distribution of particles in the plasma. The detection of 



the no n-thermal component in GX 339-4 in the LH s tate was done previously by lJoinet et al. 



(120071 ) and in the HIMS bv IPel Santo et all (l2008f ). We thus confi rm the detect i on of the 
non-thermal component in the SIMS as well (as previously done by iMalzac et al.l (120061 ) for 
Cyg X-1). 

The results obtained by applying EQPAIR fits to the data indicate a high value for 
the coronal compactness for epoch 1, but within the range of values found in the literature. 
For epochs 3, 4 and 5 this value is high as well (when compared to that obtained in epoch 
2). We thus confirm the corre lation between coro nal compactness and covering fraction of 
the cold refiecting material by iNowak et al.l (|2002| ) for epochs 2 to 5. The high value of the 
coronal compactness found for epoch 1 (HIMS) would indicate that the Comptonizing high- 
energy source iscompact in size. This would be in agreement with the proposed scenario of 
Markoff et al.l (120051 ). in which the base of the jet could be the source of the Comptonizing 
elecrons. The fact that we are detecting the thermal cut-off would be consistent with the de- 
tection of the coronal emission as well. However, for epoch 2 (SIMS), the values for both the 
coronal compactness and opacity found are extremely low. Furthermore, the kinetic temper- 
ature found for the thermal electrons in the corona is very high (and close to the high-energy 
limit indeed). We understand these as issues indicative of the absence of coronal emission 
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during this epoch. During epochs 3 to 5 (SIMS to HS), both the coronal compactness and 
opacity increase again (accompanied with the significant detection of a relativistic fine in 
epoch 3), thus indicating re-appearance of the corona after epoch 2. 

The results reported above are supported by the different fits we have done with the 
hybrid thermal Comptonization (EQPAIR) model. We fit the data with a powerlaw dis- 
tributed injection of electrons and compared with the results obtained with a mono-energetic 
distributed injection of them as well. The former distribution is expected in the case of shock 
acceleration of particles, while the second could b e achieved in reconnection events that are 
expected to power the corona (iGaleev et al.lll979l ). We found that the first model is better 
than the second for epochs 1, 3, 4 and 5. This is indicative of a corona of particles distributed 
at different speeds being the source of the high-energy emission. However, for epoch 2, al- 
beit the first proposed scenario is not discarded, and contrary to the remainder epochs, a 
mono-energetic distribution of particles is a good description of the spectrum. It seems that 
magnetic reconnection events are driving the high-energy emission in this epoch. 

Thus, we conclude that we detect spectral evolution in our data compatible with disap- 
pearance of a partis of the corona in epoch 2 (SIMS). This was followed by its re-appearance 
in epoch 3 (SIMS) and maintained in e pochs 4 and 5. Giving strength to this interpreta- 
tion is the fact that ICorbel et al.l (120071 ) detected a series of plasma ejection events during 
4-18 of February (2007) in radio, just previously to our observations of epoch 2. Also, 
the sudden increase in fiux in the 15-50 keV of the SWIFT/BAT light curve (Figure Ej) of 
15% could be related to changes occurring in the source of the high-energy emission during 
the transition from epochs 2 to 3. The possib l e disa ppearance of the corona during epoch 
2 resembles what claimed by [Rodriguez et al.l (120081 ) in the case of GRS 1915-1-105. This 
behavior could be understood as the fact that the ejected medium is the coronal material re- 
sponsible for the hard X-ray emission. The spectra of epochs 1, 3, 4 and 5 show a significant 
fraction of non-thermal particles as well, indicating that it could be due to other processes 
apart from the r mal Comptonization . For e xample, synchrotron or self- synchrotron emission 
(IMarkoff et al.l (120031 ) jMarkoff et al.l (120051 )) occurring at the base of a jet. 
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This preprint was prepared with the A AS macros v5.2. 



Table 1. INTEGRAL and XMM-Newton Observations LOG and Effective Exposure times of the INTEGRAL and 

XMM-Newton instruments. 



Epoch 


XMM-Newton ID 


INTEGRAL 


XMM-Newton 


XMM/EPIC-pn 


JEM-X 


ISGRI 


SPI 






(yyyy/ niin/dd) 


(UiC- hli:rrim : yyyy/iiiiii/ad) 










1 




2007/01/30-02/01 






50 785 


85 265 


98 324 


2 


0410581201 


2007/02/17-19 


00:03-04:44 ; 2007/02/19 


15 754 


45 036 


73 356 


85 565 


3 


0410581301 


2007/03/04-06 


11:15-11:15 ; 2007/03/05 


16 702 


49187 


93 096 


107 733 


4 




2007/03/16-18 






50 966 


82 635 




5 


0410581701 


2007/03/29-31 


14:34-20:07 ; 2007/03/30 


18 294 


63 647 


93 655 





C5 
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Table 2. Parameters obtained for the best fits of the joint XMM/EPIC-pn, JEM-X, 
ISGRI and SPI spectra (see text for details), using the descriptive phenomenological 
models constant*phabs (diskbb+gaussian+powerlaw)highecut, 
constant*phabs (diskbb+gaussian+powerlaw) , 
constant*phabs*smedge (diskbb+laor+powerlaw) , 
constant*phabs*smedge (diskbb+powerlaw) and 
constant*phabs*smedge (diskbb+gaussian+powerlaw) for the epochs 1 to 5, respectively. 

See text in Section [3TT] for more details. 



Epoch 1 



Epoch 2 



Epoch 3 



Epoch 4 



Epoch 5 



Parameter 



r 1.461 

Ecutoff ~ 

[keV] 

Ef 66. 4j 

[kcV] 

Npow 0.77 
[ph/keV/s/cm2] at 1 keV 



+0.01 



Powerlaw 
2.33±0.02 



0.89" 



0.03 



2.658^ 



5.99^ 



2.39±0.06 



0.82^ 



2.653±0.015 



0.86±0.05 



Nh 

[1022cm-2] 



0.5 (f) 



phabs 
0.467±0.008 



0.577±0.007 



0.5 (f) 



0.451±0.001 



kTin 0.97±0.03 
[keV] 

Nbb 224 
[(Rin [km] /D [lOkpc] )2 x cos0] 



+ 175 



diskbb 
0.891±0.003 
2 292^1^ 



.003 



0.852l°:°°| 0.75l:o°ol 0.698±0. 
1323j;33 3 678±1900 2 705±45 



[E] 

[kcV] 
[r] 



smedge 



;.53±0.20 



1.09^ 



1.6±0.5 



9.22±0.25 
3.0±0.5 



Egauss 

[keV] 

O" 

[keV] 

Ngauss 

[ph/cm^/s] 



gaussianj (emission ins.) 
2.28 (f) 
0(f) 
0.012±0.007 

gaussian2 (absorption ins.) 



2.25±0.03 
0(f) 
0.012±0.003 



[kcV] 
[keV] 



1.83±0.02 
0(f) 
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Tablc 2 — Continued 



Epoch 1 



Epoch 2 



Epoch 3 Epoch 4 Epoch 5 



Ng 

[ph/cm^/s] 



gauss 
.2 > 



Egauss 

[keV] 

c 
[keV] 

Ngauss 
[ph/cm^/s 



Egauss 

[keV] 

u 
[keV] 
Ngauss 
[ph/cm^/s] 



E 
[keV] 
1 

fout 

[Rg] 



N 

[ph/cm^/s 



J-^gauss 

[keV] 
<y 

[keV] 
Ngauss 
[ph/cm^/s] 



6.40+0 57 



1.83 



+0.18 



0.042 



+0.006 
-0.01 



gaussian3 (absorption ins." 



Broad gaussian4 (Fe) 



Laor 



Broad gaussiang IkeV) 
1.18±0.02 
0.09±0.01 
0.07±0.01 



6.83±0.14 



c 9+0.6 
2.11±0.11 



100 (f) 



+4 



35 



031+'' °°5 



0.010±0.004 

2.15±0.03 
0(f) 
0.007±0.002 



1.13±0.03 
0.185±0.015 
0.106±0.018 



Intrumental normalization factors 
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Table 2 — Continued 



Epoch 1 Epoch 2 Epoch 3 Epoch 4 Epoch 5 



CxMM/BPIC-pn — 
CjEM-X 1-0 (f) 

CisGRi 1.37±0.02 

Cspi 1.57±0.03 



xl 1.4 1.9 1.3 1.20 1.3 

u 91 144 141 44 102 



1.0 (f) 1.0 (f) - 1.0 (f) 

0.89±0.01 0.827±0.009 1.0 (f) 0.86±0.02 

1.39±0.02 1.18±0.02 1.28±0.02 1.38±0.02 
1.56±0.04 1.52^007 



Note. — Parameters fixed in the fits are denoted by 'f. EPIC-pn instrumental 
absorption features are implemented in XSPEC with negative gaussian profiles. 



Table 3. Luminosities (in the 0.7-300 keV energy range) and parameters (with F and being the photon index and 

the cut-off energy, respectively, see the Section [XT] for more details). 



Epoch 




Lpoui 


^total 


r 








[erg/s] 


[erg/s] 


[erg/ s] 




[keV] 




1 
2 
3 
4 
5 


5.4±1.2xl03^ 
3.1±0.7xl03^ 
1.3±0.4xl03^ 
1.3±0.2xl03s 
l.l±0.2xl03s 


2.0±0.4xl038 
4.7±1. 1x1037 
2.2±0.6xl03s 
3.4±0.4xl037 
2.6±0.5xl037 


2.6±0.6xl038 (0.33LEdd) 
3.4±0.8xl038 (0.46LEdd) 
3.5±1.0xl038 (0.45LEdd) 
1.6±0.2xl0=^s (0.20LEdd) 
1.4±0.3xl038 (O.lSLEdd) 


2.33±0.02 

2.659i°:r 
2.39±0.06 
2.65±0.02 


66.4t?:? 


0.27 
7.0 
0.6 
3.8 
4.0 



00 
o 



Table 4. Best-fit parameters of the joint ISGRI and SPI spectra for epochs 1, 2 and 3 by fitting data with 
thermal- Comptonization model and thermal-Comptonization (compps plus a powerlaw. Comptonization temperature 
(kT) and optical depth (r) are free parameters; F is the powerlaw photon-index. In COMPPS model the black body 
temperature of the seed photons was frozen to the value obtained by EQPAIR fits of XMM/EPIC-pn spectra in all 
epochs (except of epoch 2, which was frozen to 0.77 keV). Reflection component (0/27r) is only present in COMPPS. 



Epoch Model kTe t [n/27r] [r] X^(dof) F-test p 

[keV] 



1 


COMPPS 


43.6+1.4 


2.89+0.11 


0.92+0.08 




1.02 (54) 


1 


COMPPS+PO 


40.8+0.9 


2.95+0.05 


0.94+0.06 


1.1+0.6 


0.92 (52) 6.837E-02 


2 


COMPPS 


103+45 


0.48l°:?5 







1.10 (26) 


2 


COMPPS+PO 


(57+10 


0.41+0.36 





1 QQ+10 

i-yy_o.o7 


1.16 (24) 


3 


COMPPS 


135+15 


0.55+0.32 


0.37+0.25 




1.04 (26) 


3 


COMPPS+PO 


31+8 


2.95+0.05 


0.37+0.21 


1.75+0.8 


1.06 (24) 



CO 



Table 5. Best-fit parameters of the joint XMM/EPIC-pn, JEM-X, ISGRI and SPI spectra for the 5 epochs. Fits have 
been performed simultaneously with EQPAIR combined with LAOR. The inchnation angle of the disk was frozen to 
the value i — 50° for the first epoch and i — 30° for the remainder ones. See text for the parameters description. 



Epoch 








Tp 








xS(dof) 




[eV] 








[keV] 








1 
2 
3 
4 
5 


300 (f) 

704±2 
500(f) 
563±2 


n g+0.6 

0.05±Hf 
0.28in? 

Q 24+0.02 
"•^^-0.005 

r, ,0+0.012 

"•^•^-0.018 


0.90±0.10 

0.84±0.03 
49+0.02 

0.38±0.02 


2 39+0-15 
^•'-"'-0.18 

< 0.02 
1 41+0.03 
^■^^-0.06 

2.5±0.5 


27.5±1.2 

69±4 
10.8±0.3 
4.3±0.8 
10.5±0.7 


1 an+0-07 

9 07+0. 05 
^•"'-0.06 
, 42+0 03 
•^•^^-0.04 

2 Qfi+0.04 

2 64+° °® 


r, qo + 0.06 

1(f) 
1(f) 

0.72l°o;l« 


1.0(82) 
1.7(142) 
1.5(138) 
1.5(44) 
1.4(99) 
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Fig. 1. — Mosaic significance image (obtained in epoch 1) of tlie GX 339-4 region as seen 
witli ISGRI in tlie 20 — 40keV energy range, galactic coordinates from riglit to left. Besides 
GX 339-4 (detection isgnificance of 1070cr), several other high energy sources are visible. 
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Fig. 2. — Significance map obtained with SPIROS in the 200-300 keV energy range during 
epoch 1 observations; the color scale goes from -5 a to 5 o\ the grid spacing is of 1 ° and 5 ° 
in right-ascension and declination, respectively. GX 339-4 is the only source detected in the 
200-300 keV energy range, with a fiux of 0.01951±0.004030 ph/cm7sec (4.8 a). 
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Fig. 3. — Lower panels: 0.6-10 keV EPIC pn light curve of GX 339-4 for epochs 2 (upper-left 
panel), 3 (upper- right panel) and 5 (bottom panel). Upper panels: hardness ratio (counts in 
the 2-10 keV band divided by those between 0.6-2 keV) for epochs 2, 3 and 5. The binning 
is 120 s for all panels. 
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Fig. 4.— ISGRI/INTEGRAL (20-200 keV) light curve during all the observations (from 
Epoch 1 to 5). The binning is 1 000 s for all panels. 
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SWIFT/B AT (15-50 keV) 
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MJD (days) 

Fig. 5.— SWIFT/BAT and RXTE/ASM daily light curves of GX 339-4 during the overall 
outburst in 2007 (red and black dotted lines, respectively), illustrating the spectral evolution 
between the different states. Intervals of time in which the INTEGRAL observations were 
undertaken (solid black lines) and the period of time when the radio ejection events were 
detected (within both dotted black vertical lines) are also shown. 
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1 10 100 



Energy (keV) 

Fig. 6. — The plot above shows the residuals obtained when the INTEGRAL and XMM- 
Newton data of epoch 3 is fit with a phenomenological disk plus power-law model. In the 
residuals a relativistically broadened FeK^ line is present. 
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Fig. 7. — Unfolded INTEGRAL and XMM/EPIC-pn spectra (blue, black, red and green for 
the XMM/EPIC-pn, JEM-X, ISGRI and SPI data, respectively) from epochs 1 to 3 (upper- 
left to lower-right) fit with descriptive phenomenological models. The continuum line shows 
the total model (see the text Section [XT] and Table [2] for details), the dotted and dashed lines 
show the accretion disk and the powerlaw components, respectively, and the magenta line 
represents relativistically broadened FcKq, emission from a disk. The cyan colored gaussian 
lines mean instrumental XMM/EPIC-pn effects. 
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Fig. 8. — For epochs 4 and 5 spectral evolution occurred in the form of a gradual softening. 
See Section 13.11 for more details. 
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1 10 100 

Energy (keV) 

Fig. 9. — Unfolded spectra from epoch 1 to 5 (black, blue, red, cyan and green, respectively). 
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Fig. 10. — IBIS and SPI unfolded energy spectra for epochs 1, 2 and 3 (left to right) fit with 
a Comptonization model (COMPPS) plus a powerlaw. 



